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One dimensional Rectangular, Barlett, Cosine, Hann, 
Hamming, Blackman and Kaiser window are extended to two 
dimensions using circular symmetry and the Hankle transform 
in integral form. It is shown that this transform cannot 
generally be expressed in closed form. Comparisons of the 
above two dimensional windows are then made on the basis 
of two dimensional extensions of the criteria of equivalent 
noise bandwidth, coherent gain, and the so called trade-off 
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U IT LIT 


Windowing is essential to signal processing and to digi- 
tal filter design. For example, the weighting functions of 
the window are used in harmonic analysis of signals [1] and 
in design techniques for linear-phase digital FIR filter [2]. 

In fact, if an infinite sequence of numbers is truncated, 
then its Fourier transform will be subject to the well-known 
Gibb's phenomenon; which will cause over shoot and ripples 
before and after any point of discontinuity in the frequency 
domain. The overshoot and ripples will not disappear no 
matter how long the sequence is. 

To overcome this undesirable effect in the frequency do- 
main, the finite sequence should be modified using a weighting 
function, called a window. For example, if w(n) is a weight- 
ing sequence and h(n) is an infinite sequence of number, then 


the windowed finite sequence is: 
h, Cn) aa) (AO 


where: 


it OQ In|<N 


U elsewhere 


225 

w(n) = { S 
In this example the weighting functions performs both 

truncation and modification of the signal at the same time. 


Multiplication of two sequences in time domain is equivalent 





to the convolution of their transforms in the frequency 
domain. Consequently the windowing method has the effect 

of smoothing out the rough points - moderating the overshoot 
and ripples - in the original frequency response. 

The technique of windowing is shown graphically in 
ENS l.l. 

There are many window functions available in many pieces 
of literature. [1, 3] According to the way in which the 
window function is formed, they can be separated into two 
categories; basic windows and constructed windows. A con- 
structed window is formed from a product, a sum, or a con- 
volution of two or more basic windows. 

Seven best-known basic windows are listed in Table 1.1, 
and they are: Rectangular, Barlett, Cosine, Hann, Hamming, 
Blackman and Kaiser window. Most of the windows are named 
after individuals. These seven basic windows will be dis- 
cussed thoroughly in later chapters. 

Although Huang [4] shows that a good two dimensional 
window can be derived from one dimensicnal window, most 
Windows found in the literature are given in only one-di- 
mensional form. It is the purpose of this thesis to present 
and discuss two-dimensional windows derived from those one 
dimensional form, using Huang's theorem. 

The easiest way to compare various windows is to create 
some appropriate figures of merit for windows. This thesis 
presents those figures of merit that are used to compare the 


given window functions. 


LU 
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A. TECHNIQUE OF WINDOWING [2] 

A graphical interpretation of windowing is given in 
Fig. 1.1 (a non-recursive filter design using the window 
method). 

The first figure shows an infinite impulse response of 
a digital filter, h(n), and its periodic frequency response 
Bier). 

The next figure shows a rectangular window r(n), with 
its frequency response, R(w). 

he chird figure consists of multiplication of h(n) and 
r(n), which produces a truncation of the RP impulse 
response, h, (n). Multiplication in the time domain is equi- 
valent to convolution in ESTEE domain. The convolu- 
Bon of H(w) and R(w) is H Cw) containing overshoot and 
ripples in its frequency response of the corresponding digi- 
tal filter as shown in the figure. 

ME cation of impulse response is required for a non- 
recursive digital filter design, since the filter produces 
an output which depends only upon the present and a finite 
number of past inputs. 

The fourth figure shows the non-rectangular window w(n), 
and its periodic Fourier Transform W(w). The purpose of 
this window is to alter the truncated filters impulse response, 
h, (n), slightly in order to minimize the overshoot and ripples 
in its frequency response H (w). 

The last figure shows the finite impulse response h(n), 


Beh consists of multiplication of h, (n) and w(n). The 


ih 





frequency response of h,,(n) is produced by convolving H. (w) 
with W(w), resulting in a smooth frequency response H (w). 
The technique of windowing applied to a filter design, 


results in an improved filter with a smooth frequency response. 


B. THESIS ORGANIZATION 

The following discussion presents the organization of 
this thesis. 

As mentioned before, the best way to discuss windows is 
to define parameters which can be used to characterize a 
window.  Logically, the window figures of merit are first 
discussed in the next chapter, Chapter II. 

Chapter III gives an overview of one-dimensional windows 
and specifies their characteristics, using parameters defined 
Em Chapter II. Table IIl.l, given at the end of this chapter, 
tabulates one-dimensional windows and their parameters. 

The transformation from one-dimensional to two-dimensional 
is given in Chapter IV, where we present three major discus- 
sions. First we present the expansion of the Fourier Trans- 
form from one-dimension to two-dimensions. Then we present 
how Huang [5] expanded a good one-dimensional window into a 
good two-dimensional window. Finally, the definition of one- 
dimensional figures of merit are expanded into two-dimensional 
window figures of merit. 

Chapter V presents two-dimensional windows in detail, and 
specifies their parameters as defined in Chapter IV. Table 
V.1, given at the end of this chapter, tabulates two-dimen- 


sional windows and their parameters. 
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Figure l.l 
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Application of Window Function To 


FIR Filter Design 


de 


IT 





Two-dimensional Fourier Transform of window functions 
are hard to present in a closed form. Fortunately a closed 
form Fourier Transform is not required in studying windows, 
so two-dimensional Fourier Transform of the windows are 
presented in integral form. The effort to present a two- 
dimensional Fourier Transform of a window in a closed form 


is given in the Appendix. 
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RE 


Seven Basic One-Dimensional Window Functions 


No. Window Punction 
1 Rectangular W(x) = 1.0 |x| <a 
2 Barlett ico gr nr |x| <a 
3 Cosine W(x) = cos T x Ix| «a 
4 Hann W(x) = H : 5 cos = E |x| <a 
5 Hamming W(x) = 0.54 + 0.46 cos X |x| <a 
6 Blackman W(x) = 0.42+0.5 cos 3 x + 0.08 cos 
|x| <a 


AE 
113 d 121020) } 


ET le 
Ig B 


7 Kaiser W(x) 


jo 


pa 





I TRT [TEIR [1, 3] 


Since windows are used in frequency response analysis to 
reduce undesirable effects - overshoots and ripples in the 
frequency domain - it is very important to find the para- 
meters or figures of merit of the windows that can be used 
to identify their characteristics and to indicate their per- 
formance. 

There are many convenient measures that can be used as 
Parameters. 

The common desirable characteristics of a window are: 

l. A narrow main lobe width of the frequency response 
of the window, containing as much of the total energy as 
possible. 

2. Side lobes of the frequency response decreasing in 
energy rapidly as frequency increases. 

These two common characteristics certainly can be used as 
parameters to identify the window's performance. 

Harris [1] and Geckinli [3] define one-dimensional window 
performance parameters which will be repeated here. These 
parameters are used throughout this thesis and will be expan- 
ded for the two-dimensional case in later chapters. The 
parameters are: equivalent noise BW, coherent gain, and 


trade-off parameters. 
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A. EQUIVALENT NOISE BANDWIDTH (ENBW) 

The ENBW parameter is the width of a rectangular filter 
in the frequency domain with the same amplitude as the peak 
power gain of the window, and with a width such that the 
total noise power in the rectangular filter is equal to the 
total window noise power. 

A graphical explanation of the parameter is shown in 
Fig. 2.1. The first figure shows a window frequency response 
which consists of a main lobe and side lobes. The second 
figure shows a discrete frequency component of a signal, 
together with white noise. The third figure shows the result 
of windowing. 

The window's frequency response convolves with the dis- 
crete component of the signal plus noise. In this case the 
window acts as a filter, gathering information over its 
bandwidth. The last figure represents the area of the windowed 
discrete frequency response of the signal plus noise, at dis- 
crete frequency ue 

It is clear from the figure that to minimize the accumu- 
lated noise collected by the window, it is required that the 
window has a narrow bandwidth. 

A convenient measure of window's bandwidth is called 
equivalent noise bandwidth. 


Noise power accumulated by the window is defined as: 


oo 


Noise Power = — J. ww) |? dw C27 A) 





where Jc noise power per unit BW. 
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Spectral Amplitude 







noise 


MIME MA k 





Figure 2.1 


Graphical Interpretation Of Convolution Of A Window 


And Spectral Amplitude Of A Signal, At ww 


TE 





Parseval's theorem allows equation (2.1) to be computed as: 


N T 2 
Horse rovar -= m "EU o» dt (252) 
-T 


The integral in (2.1) and (2.2) represents the area under 
the window curve. If the area is divided by height, then 
the result is the width of a rectangle that has the same area 


Dh the window. 


The peak gain of signal occurs at w=0, and is defined as: 


Hr 
Peak signal gain = w(0) = / w(t) dt RJ 
-T 
: 2 A 2 
DT = ow (oO) =[ f w(t) dt] (2.4) 
-T 


Then the ENBW is equal to the noise-power divided by peak 


power: 


I Z 
ee) dt 


~ -T 
MEM alo 3 (2,5) 


RR a dt] 
ET 


The above expression is a normalized ENBW, after it is 


divided by factor o . This parameter will be used as a 
2T 


window figure of merit. 


Bee COHERENT GAIN 

Coherent gain is best defined using a sinusoidal signal. 
In Fig. 2.1 a component of signal frequency response at a 
discrete frequency is shown. Convolving this signal with a 


window's frequency response can be thought as passing the 


103 








Peak Power Gain = Iwco)]? 


Equivalent Noise BW 
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Figure 2.2 


Equivalent Noise BW Of A Window 


20 





signal through a filter. Using this view, the coherent gain 
of a window can be examined by passing a sinusoidal signal 
through a filter. 


Let the sinusoidal signal be defined as: 


mer) = AeJYX TEC) (20959 


where n(x) - white noise. 


The filter output on the windowed signal is given as: 


f(x) = A Ere T + n(x) w(x) - a<x<a 7) 


The frequency response of the windowed signal only (less noise) 


is given as: 


a : : 
F (£F) = J A W(x)eJ"* g7J"* ae = 
W 
= 
a 
= A f W(x) dx (2.8) 


=a 


The windowed frequency response is proportional to the input 
amplitude A. The proportionality factor is the integral of 
the window function. 

For a rectangular window this factor is 2a. For any other 
window the factor is reduced due to the window smoothly going 
to zero near the boundaries. 

The coherent gain is defined as the proportinality factor 
divided by factor 2a, the extension of the window (normalized). 
a 


f W(x) dx (2.9) 
-a 


COH.GAIN = — 
2a 


Zl 





The coherent gain is the second parameter that will be used 


as window figure of merit. 


C. TRADE-OFF PARAMETERS 

Trade-off parameters are best explained using Fig. 2.3. 
The top figure shows the Hamming window in the time or spa- 
tial domain, and the bottom figure shows its Fourier trans- 
form, which is given in log frequency versus normalized log 
magnitude (20 log W(f)/W(0)). 

The window's characteristics in frequency domain are 
largely determined using these trade-off parameters. 

From Fig. 2.3(b) these parameters are: 

E. aj: highest side lobe level in dB. 


S a, i side lobe level at a higher frequency, 
where asymptotic assumption is valid for 
the window behavior. 

Ser D : frequency at which main lobe drops to the 
peak side lobe. 
4. d : side lobe fall-off rate in dB/octave. 
Clearly from the above definitions, the smaller the value 
for b, the better the window performance; the smaller the 
value for aj» the smaller the leakage through the near side 


lobes; the smaller the value for a and the larger the value 


2 
of d, the smaller the leadage through the far side lobes. 
These parameters are expanded for two-dimensional case in 


a later chapter. 
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Plot Of Hamming Window, Is Used To Illustrate The 


Definition Of Trade-Off Parameters a1» an Brand d. 
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III. ONE DIMENTIONAL WINDOW FUNCTIONS 
A 7] 


The seven common window functions listed in Table I.l, 
are reviewed in detail in this chapter. Figures of merit 
for each window are presented using the definitions given in 
the preceding chapter. The window functions treated here are 
in continuous form rather than in discrete form. 

Without loss of generality, let the window function, 
w(x), be unity at the origin and limited to the interval 


|x|< a. Presented in mathematical form this means: 


PRM Cr) di 2-1 EE 


-00 


w(o) 


w(x) = 0 for |x|>a (252) 


The Fourier Transform of the above window is given as: 
a 


e ee, |x| <a CESEN 
-a 


The transforms of the windows generally consists of a main 
lobe with a large peak amplitude and side lobe with compara- 
tively small peak amplitudes. To better display the presence 
of the side lobes, logarithmis scales are used for plotting 


ın place of linear scales. 
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A. RECTANGULAR WINDOW [1] 

The rectangular window is commonly called a boxcar or 
Dirichlet window. Its main purpose is to truncate the 
length of a function to a finite extent, without affecting 
its amplitude. 


The expression for this window is given as: 


w(x) = 1.0 |x| <a (3.4) 


a 
e 


The Fourier Transform of the window is given as: 


Wf) - 3 f e"* ax 
a 
-a 
Sin3W 2 sinc(aw) (3.5) 
aw 
Ns rPansform is seen to be a Dirichlet kernel. In this 


report the extent of the window is taken as a=0.5, and the 


ENBW is given as: 


05,5 
H dx 
= 720.5 e 
ENBW = No 5 = 7120 
[ fs dx ] 
-0.5 


The coherent gain is given as: 
Ds 


COH.GAIN = f day «bs 
-0.5 
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Plots of the window function and its tranform are given 
Eric. 3.1(a) and (b) respectively. The main lobe width 
(the width between origin and the first zero crossing) is 
1.0, and the first side lobe peak is approximately 13 dB 
down from the main lobe peak. The side lobe fall off is 


6.0 dB/octave. 


B. BARTLETT WINDOW 

The Bartlett window is also called Fejer or triangular 
window, and can be thought of as a convolution of rectangu- 
lar window with itself. This window is the simpliest of all 
Windows that changes the data values by multiplication by a 
number other than unity. The mathematical expression of this 


window is given as: 


w(x) = 2.0 - X |x| <a EN) 


a 
The transform of this window is given as: 


e S 
WC£) = = O Ll ¿du dx = 
= G 


O 





. E : 
E (1.0 = “je J"* ax + = EP (1.0 = 20 I™* ax = 
a a a a 
-a O 
EU? -iux ab > -IWX 1 8 -3wx 
ER f e J dx + za f xe J dx t > fe J dx - 
EI aà 
a a` -a O 
a Em 
] fx e J dx = 
2a O 
BE. -jwx O 1 9 ENS 1 ,e7JWX a 
ES -jw PEE, x jiw ag pa 
-åa a -a 
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FReure 3. 


(b) Spectral Window 


(a) Rectangular Window, 
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J de 
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a a JW a2 la 
a -Jjwx 
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E. > E +3aw ENSURE M -1 (e JAN o + oe Jaw, 
aw a w 
2 jaw/2 _-jaw/2 - 
ME) es, 
aw/ 2 23 
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EE 27 - (ine av/2)' SEH 


The transform is seen to be the squared of Dirichlet Kernel. 
If the extent of the window is taken as a-0.5, then the 


parameters ENBW and COH.GAIN are: 





.-0.5 MEUS - 
ENBW z NOMEN c 7? = Se je 
[ sf (1- = dx J] 
EN (DES 
0.5 = 
COH.GAIN = ff CL =- ——) dx = 0.5 
-0.5 WE 


The window function and its transform are shown in Fig. 
3.2(a) and (b) respectively. Its main lobe width is 2.0, 
twice that of the rectangular window, and the side lobe fall- 


Off is 12 dB/octave. 


fee COSINE WINDOW 
Cosine window is actually one member of a family window 


with a general form given as: 


w(x) = cos” (ex) |x| <a (3.8) 


Type of this window depends upon parameter a, which usually 
has a value between 1 and 4. Of particular interest in this 
family window is Hann window, where a=2, which will be con- 
Sidered later. 


The cosine window is defined with a=l, as: 
w(x) = cos zx) |x| <a (3.9) 


The Fourier Transform of this window is given as: 


T - wx 
cos —x e J dx 


w(f) = 23 


D Za D 


ES 
2a 
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Figure 3.? 


(a) Barlett Window, (b) Spectral Window 





a ! T = T 
E Jt Pal CUR E dx * ES f e) Cut) x dx = 
2a E E- 


. TT . TT 
= c -— — —À i 
Sin a(w-55) sin alw ja) 
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n5| — 
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Ne cransform of this window is a summation of two shifted 
Dirichlet Kernel. 
If the extend of the window is taken as a-0.5, then the 


parameters ENBW and COH.GAIN are given as: 





(S 2 
y (costx) dx 
NES es 

— - 055 NEUES 2 

[ f costx dx] 

-0.5 

0.5 2 
eoD. GAIN = Sf cosıx dx = = 0.64 
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Plots of the window and its transform are given in Fig. 
3.3(a) and (b) respectively. From the transform plot, the 
main lobe width is 1.5 times that of the rectangular window. 
Ihe first side lobe peak is approximately 23 dB down from 
the main lobe peak, while the side lobe fall off rate is 


IENdgB/octave. 


e. 
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D. HANN WINDOW 


The Hann window is one member of the family windows given 


by cos x; where a=2. The window equation is written as: 
n eme. T a TX 
HU» cos == 7 t cos — |x| <a (SAL) 


The Fourier Transform of this window is given as: 


E à 
RCE) = = Pos as Bel) am dx = 
E al 
m 
a : d f 
Em 1 -J)WX l T =JwWX E 
^ mm E e dx + TE Sé mos EU e Goad. = 


sin a(w-L) sin a(2+1) 
a a 


l sin aw 
2 


+ À 
. aw 4 


+ 
Fir 


= C RT?) 
a(w--) 


a(w*L) 

a 
The Transform of this window consists of a summation of 
three Dirichlet Kernels, one at the origin and the two others 
shifted = radians from the origin. 

The three terms on the transform of this window are 
shown in Fig. 3.4, where the shifted Dirichlet Kernels are 
located at the first zeroes of the center kernel. The peak 
amplitudes of the shifted kernels are half that of the center 
kernel; also the side lobes of the translated kernels are 
about half the size of, and of opposite phase from those of 
the center kernel. The summation of the three side lobes in 


phase opposition, tends to cancel the overall side lobe structure. 
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Figure 3.4 


Transform Of Hann Window As A Sum 
Of Three Dirichlet Kernel 
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The cancellation of the side lobes structure suggests 
a constructive technique to define new windows. The best- 
known windows with this type of side lobes cancellation are 
Hamming and Blackman windows, which are given in the next 
section. 

If the extent of the Hann window is taken as a=0.5, the 


parameters ENBW and COH.GAIN are given by: 


o rs cos)? dx 
_ -0.5 E Ole hs 
ENBW = A a o > RES 
L f (3t+% cos Ex) dx] 
-0.5 i 
0:5 GER 
ESH-GAIN = Sech (45445 cos om: DT 0.5 


The window and its Fourier Transform are plotted in Fig. 
3.5. The main lobe width is 2.0, twice that of the rectan- 
gular window and the first side lobe level is 23 dB down 
from the main lobe. The side lobe level fall off rate is 


18 dB/octave. 
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Figure 3.5 


(b) Spectral Window 


(a) Hann Window, 
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E. HAMMING WINDOW [5] 

Hamming window can be thought an an improvement of Hann 
Window, in a sense that the first side lobe cancellation 
is even greater. From Fig. 3.4, the summation of the three 
side lobes produce an inexact cancellation. The cancellation 
can be made more pronounced by adjusting relative size of 
eaeh term. 

Exact cancellation of the side lobes was first shown by 


Hamming using a window given in general form as: 
w(x) = a + (1-0) cos = x |x| <a Cos) 


ee is a coefficient that will be defined later. The 


Fourier transform of the. window is given as: 


meee) = asinc aw * Q.5(l-a)sinc a(w-5) 


HOL Skla) sinc a (w+) (3.14) 


EE (03.147) the best cancellation of the first side lobe 


@@eurs when a = SE : = = 0.543478261). If a is selected as 
0.54, an approximation of ZE then an improvement in cancel- 
lation of the side lobes occurs. The value a=0.54, gives 


rise to a window called the Hamming window, which is written 


as, 


w(x) = 0.54 + 0.46 cos 7 . |x| <a (3.15) 
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The Fourier transform then is given as: 
Wit) = 0.54 sinc aw + 0.23 sinc a(w-—) + 
LI sme a(w+-) À (3716) 


If the extent of the window is a = 0.5, then the para- 


meters ENBW and COH.GAIN are given as: 


o s T 2 
J (0.54+0.46c0s EX) dx 
_ -0.5 i 0.3974 
ME Gs Date "7 
Es (0.54+0.46c057X) dx] 
-0.5 : 
OS - 2 
T (0.54+0.46cos—.X) dx 
EEE ` T as 
COH.GAIN = ME ee ot - 0.2916 ` SÓ 
L 7S (0.5440 .46cos5——x) dx] 


The Hamming window and its transform are ploted in Fig. 
3.6(a) and (b) respectively. From its transform plot, the 
main lobe width is twice that of the rectangular window, and 
the first side lobe level is 43 dB down from the main lobe. 


The side lobe fall of rate is 6 dB/octave. 


F. BLACKMAN WINDOW [6] 
Hann and Hamming windows are examples of windows con- 
Ene ted as the summation of shifted sinc(x) functions. 
Blackman suggested a general rule by which window 


Ame tions could be constructed as follows: 
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Figure 3.6 


(b) Spectral Window 


(a) Hamming Window, 
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a, Cos m (=) |x| <a (3.17) 
and subject to the constraint: 

a - 1.0 (3.18) 
The Fourier transform of this window becomes: 


k m a 
W(f) - 2 (-1) 


= {sine a(w-—m) + sinc a(wtom)} Cored} 
m-0 3 = 


Hann and Hamming windows are of this form with a and a 


UN es 


being non zero. Their spectral dw are summations of 
three shifted kernels. 

Any value of k will produce a new window with a summation 
of (2k-1) kernels in its tranform, and k non zero corficients, 
however one way to A windows with a narrow main lobe 
is to restrict k to a small interger. 

Blackman examined this window with k-3, and found best 


cancellation of the third and fourth side lobes occurs when: 











8j 7 iggog ^ 0:525 590 71 8 0.42 
NE a | 

a, = mg ^ 0.496 560 62 = 0.50 
Bo. s 

a, = vgggg ^ 0.076 848 67 3 0.08 
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IUT T T whieh uses these three coeficients 


is known as the Blackman window, and is given as: 


w(x) = 0.42 + 0.5 cos + 0.08 cos—x |x| <a en) 


ne Fourier transform of this window is: 
PRESS = 0.42 sinc aw + 0.25 eine a (w-5) Cs cine a(wtt) 


RG O! cine a(w-21)+0.0u sinc a (ie AT) 
Co 210) 


If the extent of this window is a-0.5, then the ENBW and 


COH.GAIN parameters are given by: 


[4.5 = 27 2 
ORO 5 cos EU * 0.08 cos EU dx 
ENBW = o 5 > E 
Ree Ss cosix + 0.08 cos EET dx] 
ns = 2 
~ 0.2646" 
Se D 
SE T Zu 
COH.GAIN - ed EE cos 45 + 0.08 cos E dx pa0e 06 
-0.5 


The window and its transform are shown in Fig. 3.7(a) and 
(b) respectively. The main lobe width is 3.0, three times that 
of the rectangular window, and its first side lobe level is 
41 dB down from the main lobe. The side lobe fall off rate is 


EwdB/octave. 
41 
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Figure 3.7 


(b) Spectral Window 


(a) Blackman Window, 
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G. KAISER WINDOW [7] 
The Kaiser window is also called a Kaiser-Bessel window 


and its mathematical expression is given as: 


I [8 /1.0-(&)°3 


a(k) = EE RS |x| «a (222) 
O 
where I (x) = zeroth order modified Bessel-function, and is 


@eeined as: 





ieee x) = 1 + ] Can Ze) 
O 


nm 8 
r1 


The Fourier transform of this window is given as: 


SE Buc 


O EZ 
B n se 


W(£) f <fa CSch) 


where fa = 8/mT. 


From (3.24) if f=fa, the term inside the square root 


becomes zero and W(fa) = E If f>fa, the term becomes 


U 
) £ 2 | EMEND. 
iy (75) ameno he identity sinh ix = à sin x, (3.24) 
becomes: 


sin[ß v (2) 2-1] 


w(f) = EE MM S M llo AR Dd CASO, 


CR 
2 SÉ FDA 
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The usual range of 8 is between 4 and 9, corresponding to 
a range of relative peak heights from 3.1% down to 0.047%. 
As 8 increases, the side lobes decrease in amplitude while 
the main lobe widths increase. By varying 8 the proportion- 
ality between the width of the main lobe and the height of 
the side lobe can be adjusted as desired. 

When it is given ß=6 and a=0.5 then the parameters ENBW 


and COH.GAIN are given as: 


Se eA E E E 


Bye = | 
ENBW = ais A -. 0.3667 = 1.5 


0825028 
S r 6/1- Bm 


L b. — and 


ES 
0.5 E 


COH.GAIN - J —— dx = 0.5 
= 09 1,16] 


Fig. 3.8(a) and (b) shows the window and its transform 
respectively. The main lobe width is 2.0, twice that of 
the rectangular window, and the first side lobe level is 
43 dB down from the main lobe. The side lobe fall off rate 
is about 6 dB/octave. 

For B=6 and 8=9 the window resembles to Hamming and 


Blackman window respectively. 
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Figure 3.8 


(b) Spectral Window 


(a) Kaiser Window, 
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ONE-DIMENSIONAL WINDOW AND FIGURES OF MERIT 
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IT TITI TT IER ONE-DIMENS ION 
TO TWO-DIMENSION (4, 8, 9] 


The two-dimensional Fourier transform, two-dimensional 
window functions, and window figures of merit are discussed 


in this chapter. 


A. TWO-DIMENSIONAL FOURIER TRANSFORM [8] 

The Fourier transform of f(x,y), a function of two inde- 
pendent variables x and y, is represented by F .{f(x,y)} or 
F(u,v), as is defined as: 


co oo 


Bey} = FOnv) s f£ / fu ye UE UY daddy —— (1) 


== 00 = OO 


The transform is itself a function of two independent 
variables u and v, which are generally referred to as fre- 
quencies. 

Similarly the inverse Fourier transform of a function 
F(u,v) is represented by F + {F(u,v)} or f(x,y) and is 


defined as: 


rov) ray), ip dom a 2) 


As mathematical operations, the transform and its inverse are 


very similar differing only in the sign of the exponent 
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appearing in the integral. 
In order that the integral exist, the constraints on 
Eme function f(x,y) are given as: 
l. f(x,y) must be absolutely integrable over the infinite 
Xoy plane. 

2. f(x,y) must have only a finite numbers of discon- 
tinuities and a finite number of maxima and minima 
in any finite extent. 

3. f(x,y) must have no infinite discontinuites. 

The basic definition of the Fourier transform pair leads 
to a rich mathematical structure associated with the trans- 
form operation. Considered below, a few of the basic mathe- 
matical properties that will find use in later discussion. 
These properties are presented as mathematical theorens, 
followed by a brief statement of their physical significance. 


l. Linearity theorem: 
E E F 
Eene B g) (4.3) 


The transform of a sum of two functions is simply the sum 


EE cheir individual transforms. 


2. Sealing theorem: 
PCa Diy) | = gr G (=, ei (4.4) 


A stretching of the coordinates in the space domain (x,y) 


Meemicts in a construction of coordinates in the frequency 
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domain (u,v), plus a change in the overall amplitude of the 


speetrum. 
EE Shifting theorem: 


F {f(x-a, y-b)} = Bea ca? Co) 


A translation of a function in the space domain intro- 
duces a linear phase-shift in the frequency domain. 


4. Parseval's theorem: 


GO GO co 00 


D "ELE my : dxdy = f J |F(u,v)|* dudv (4.6) 


=D -CO - COO = CO 


This theorem is interpretable as a statement of conserva- 


tion of energy. 


5. Convolution theorem: 
pei J ED -xdezdnb s FCu,v) G(u,v) (ue 


The convolution of two functions in the space domain is 
entirely equivalent to the simpler operation of multiplying 
their individual transforms. 


6. Autocorrelation theorem: 


ae. f if BE (er y n)dzdn) ee (4.8) 


49 





similarly: 


CO 


r {|f(x,y) |f} E J F(z,n)F (z*u,n*v)dudv (4.9) 

These theorems are special cases of the convolution 
theorem. 

fe Fourier identity: 

Peach point of continuity of f(x,y): 


Ft ve (f£(x,y)) * F EF 3 (f(x,y)0) - £(x,y) (4.10) 


The successfive transformation and inverse transforma- 
EU of a function yields that function again, except at 


Peame Of discontinuity. 


1. Separable and Circularly Symmetric Functions [9] 


A function of two independent variable is called 
separable with respect to a specific coordinate system if it 
can be written as a product of two functions. 

A function f(x,y) is separable in the rectangular 


coordinate x,y if: 
imo ,y) = f(x) f, (y) (4.11) 


A function f(»r,8) is separable in the polar coordinate 


foe if: 


f(r,8) = f,(r) f,(9) (4.12) 
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Separable functions are more convenient to deal with 
than more general functions. Separability allows complicated 
two dimensional manipulations to be reduced to simpler one 
dimensional manipulations. 

A function separable in rectangular coordinates has 
the simple property that its Fourier transform is obtained as 
a product of two one dimensional transforms 


GO GO 


RE a 


Peer (x,y) } 


f Bern e: f f(yle dy 


«= 00 


F (É, (x)) F Way ES 


Thus the transform of f(x,y) is itself separable into 
a product of two factors. The process simplifies to a familiar 
one dimension manipulation with a function of u, and a function 
SUV. 

Functions separable in polar coordinates are not so 
easily handled as those separable in rectangular coordinates. 
But it is generally possible to demonstrate that two dimen- 
sional manipulations can be performed by means of a series of 
one dimensional manipulations. The simplest class of functions 
separable in polar coordinates are those functions which ex- 


hibit circular symetry. 
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A function is said to be a circularly symmetric if it 
can be written as a function of. radius r, independent of the 


Mariable 9. 
He. 0) = f(r) aa, 


Such functions are very important, since a window function can 
be extended into a two dimensions 1f it is given a circularly 
symmetric form. Accordingly, special attention should be 
Pema to these functions. 

The Fourier transform of a function in rectangular 
coordinates is given by: 


00 oO 


BEE ty) xay Gis) 


« 00 we OO 


If both x,y and u,v coordinates are transformed into 


a polar coordinate system then: 


r = L A x=rcosd 

GE = tan" LÉI y=rsin9 
and 

N: Au u-pcosóo 

9 = tant) v=psing 


Substituting the polar coordinates into (4.15) and 


using the transform identity: 


SI 





ady = pagar 
uxtvy = pcosgr cosbtpsindr sind = r cos(9-$) 


emos er eransiorm in (14,15) becomes: 


p zu -jor cos(8-9$) 
Pio Ais P Ss rn f(r,e)e 1? " 48 dy (4.16) 
U U 


If the Bessel function identity: 


2T : 
ERC or) = ES mes Cs cos (6-9) dé Cue) 
O 2T 0 


is used to simplify (4.16), then the transform becomes: 


CO 


F(p) - 2m d PE (r) J Cer) dr (4.18) 
Thus the Fourier transform of a circularly symmetric 
function is itself circularly symmetric, and can be found by 
performing a one dimensional manipulation as per (4.18). The 
expression in (4.18) is known as the Fourier-Bessel transform 
or as the Hankel transform. 
By identical arguments, the inverse Fourier transform 
of a circularly symmetric function F(p) can be expressed as: 


CO 


mer) = 27 $ pE(p) J (rp) dp (919) 
U 


Se 





There is no difference between the transform and its inverse 
operation in the case of Hankel transform. Using the notation 
H{ } to represent the Hankel transform operation: 


MEE) E EA Eto) = fir) (4.20) 


From the Hankel transform pair (4.18) and (4.19), two 
dimensional Fourier transform may be found from one dimen- 
sional mainpulations. 

Huang [4] has shown that a good two dimensional window 
function can be derived from a good one dimensional function, 
as long as it is circularly symmetric. This is discussed in 
the next section. 

The Hankel transform in (4.18) will be used to cal- 


Ate the Fourier transform of a two dimensional window. 


B. TWO DIMENSIONAL WINDOW FUNCTION [4] 
Ef W, C(u,v) and Wo (x,y) is a two dimensional Fourier 
transform pair, then its inverse transform relation is given 


as: 


oo 00 


we (x,y) = f Jr W Qu, v)el C * VY) avd (4.21) 


For y=0, (4.21) becomes: 


co o 


Ww, (x0) ee W,(u,v) dv}e 


- CO CO 


JUX du (108225 


ol 





TE w„(X,y) is given as a circularly symmetric function, 


then Wo (XY) = Tee Perper reularly symmetric function 
w,(x,0) BUCO Ci 23) 


And its inverse Fourier transform is given as: 


me) = ff Wiel 


a= CO 


du C24) 
From (4.22) and (4.23) it is evident that: 


Wu) = ff W,Cu,v)dv (4,25) 


=% 


The expression in (4.25) is called a profile function, 


and similarly the profile function in v is expressed as: 


CO 


W(v) - f W,Cu,v)du (4.26) 


een CO 


Now consider the following two equations: 


E or O 
pou) - íi Cu a) 
U Wc o 


and 


i neu o for all v 
H, (u,v) x (127195 
U E all v 
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The expressions in (4.27) and (4,28) are one-dimensional and 
two dimensional step function respectively. 

LE W, (u,v) is a two dimensional window's Fourier trans- 
form, then its convolution with the step function is given 


a5. 


W,(u,v)8 H, (u,v) = f D W,(t,s)H(u-t, v-s)dtds 
u 00 
er X Wi (t,s)ds}dt (4.29) 
O — 00 


where t,s are dummy variables. 
Using the profile-function identity given in (4.25) and 


REDE), it can be shown that: 


CO 


W(t) = / W,(t,s)ds Goen 


== 00 


Equation (4.29) becomes: 
u 
W,(u,v) 9 H, (u,v) ew tat (4.31) 
O 


Bepending the right hand side of (4.31) into a convolution form: 


CO 


u 
PCE 3 dT f W(t) H(u-t)dt 
O es) 


Wa) 8 H(u) (4,32) 
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Nu substituting (4.32) to (4.31): 


W, (u,v) 9 H, (u,v) = W(u) 8 H(u) een 


Huang states that if a two dimensional function is given 
in circularly symmetric form, then its Fourier transform 
satisfies the relation given in (4.33). 

To show the benefit of the results in (4.33) consider an 
example of designing a low pass filter. 

EXAMPLE : 
The ideal Fourier transform of one dimensional Low Pass 
Filter is given as: 
1 for jul < A 
meu) = i Cre) 
0 jul > A 

The impulse response f(x) of the filter is equal to its 
inverse Fourier transform, which has an infinite duration. 
In reality, the impulse response has to be finite. By using 
a window function, the infinite duration of the impulse res- 


ponse is truncated to become a finite impulse response as 


given as: 

g(x) = f(x) - w(x) (4,35) 
where: 

ax) = 0 mon! oda. SEN 


oy 





The windowed frequency response of the FIR low pass filter 


is now equal to the convolution F(u) and W(u): 
eu) = F(u) 9 W(u) Qucm) 


The ideal Fourier transform of a two dimensional LPF is 


given as: 


1 |u^*v? | «à? 
0 Ae 


Fu, v) = { (4.38) 


The infinite impulse response f,(x,y) © this filter has 
to be truncated using a window function w„(x,y). The finite 


impulse response filter is given as: 
g, (x,y) = f(x,y) wy (x,y) Gis os) 
The windowed frequency response of this filter is given as: 
G, (u,v) = F C(u,v) 8 W, (u,v) (4.40) 
A good window should have the width of its frequency response 
much smaller than the width of the filter's frequency response. 


W(u) and W,tu,v) should be zero for values of Kees much 


smaller than A. 
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Near a discontinuity in the filter's frequency response 
Be 1.00, u = tA, Glu) 1s essentially a convolution of 


W(u) and a step function: 
G(u) = H(u) 8 W(u) (4,41) 


Similarly near a discontinuity in the filter's frequency 
response F (u,v), d. a = G,(u,v) is essentially a 


convolution of W,(u,v) and a step function: 
G, (u,v) = H, (u,v) 9 W, (u,v) (4.42) 


From Huang's theorm given in (4.33), G,(u,v)=6(u) if only 
we (x,y) = w(y x a 

This means that if w(x) is a good window, then, w„(x,y) 
is also a good two dimensional window. 

If other words, a good two dimensional window can be 
derived from a good one dimensional window by extending 


BEI a circularly symmetric function. 


C. TWO DIMENSIONAL WINDOW FIGURE'S OF MERIT 

The window's figures of merit found in the one dimen- 
sional case are now applied to the two dimensional case. The 
same figures of merit that may be used for comparison among 
Windows, will be re-examined to determine if there should 


be any changes for the two dimensional case. 
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1. Equivalent Noise Band Width (ENBW) 

The ENBW in one dimensional case is defined as the 
width of a rectangular filter with the same peak power gain 
that passes the same noise power. 

A rectangular filter in two-dimension represents a 
cross-section of a cylindrical filter. The noise power 
passed by the area under the rectangular filter, becomes 
noise power passed by the volume under the cylindrical filter. 

If the peak power gain is the same, then the volume 
divided by the height, results in the area covered by the 
base of the cylinder, which is a circle. 

Two-dimensional ENBW is defined as the circumference 
of the horizontal cross-section of a cylindrical filter, 
that passes the.same noise power. 

The vertical cross-section of a two dimensional 
window is equal to a one dimensional window. Consequently 
the width of the equivalent rectangular window is equal to 
the diameter of the equivalent cylindrical window. Then 


f 


the two dimensional ENBW is given as: 
(ENBW) ^p = T(ENBW)- p (4.42) 


For the two dimensional case, the definition of ENBW 
should be changed from the width of a rectangular to the 
circumference of a cylinder, and the value can be calculated 


directly from the one dimensional case. 
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2. Coherent Gain 
Two dimensional coherent gain is best defined using 
two dimensional sinusoidal signal. The input signal is 


defined as: 


SI) (4.43) 


hee) = me 


where n(x,y) is two dimensional white noise. 


The output signal or the windowed signal is given as: 


3Cux*vy) 


f(x,y) - Aw(x,y)e tw(x,y)n(x,y) (4.44) 


The two dimensional frequency response of the win- 


dowed signal is given as: 


Ee l 
HEC) fo JT pr e eee a, 
ii -a -a 
t cl 
=A S DTE EON. (4.45) 
-a <a 


The amplitude of the windowed signal frequency response 
is proportional to the original amplitude A. The propor- 
tionality factor is the double integral of the window function. 

If the window function is given in circularly symme- 


tric form, then the proportionality factor becomes: 


a a am 27 
J ana y ddy = J Aaa 
-a -a O O 
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a 2T 
A ON edr J do 


O O 
a 

- 27 frwl(r)dr (4.46) 
O 


The proportionality factor for a rectangular window 
is equal to ra”. The proportionality factor for the other 
window is reduced due to the window smoothly going to zero 
near the boundaries. 

If the proportionality factor in (4.46) is normalized 
by factor ma’ (that 1s, divided by the area covered by the 
base of the window), then the proportionality factor becomes 
coherent gain: 


BOHSGAIN = pr w(r)dr GE) 


O == py 


e 
2 
a 


The coherent gain given in (4.47) above is used as 
window's figure of merit. 
3. Trade-off Parameters 
Two dimensional trade-off parameters are similar in 
definition with the ones used in the one dimensional case. 


The trade-off parameter a b and d are derived from the 


Br 
frequency response of a window, which is found using the 
Hankel transform. 

The Hankel transform allows a two dimensional trans- 


form to be found from a one dimensional manipulation.  Con- 


sequently the definitions of the trade-off parameters doesn't 


62 





nge. The graphical interpretation of the parameters are 


Mewn in Fig. 4.l. 
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Figure 4.1 


Plot Of Kaiser Window, Used To Illustrate The 


Definition Of Trade-Off Parameters aj» 855 b and d 
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V. TWO DIMENSIONAL WINDOW FUNCTIONS 


The two dimensional extension of the previous window 
functions will be given in this chapter. The window's fre- 
quency response is found using Hankel transform.  Unfor- 
tunately it is hard to present the Hankel transform of a 
window in a closed form. Only two transforms are presented 
in closed form, they are Rectangular and Barlett windows, 
the rest are presented in their original integral form. 

Without loss of generality, the window function is given 
as w(r), a two dimensional circularly symmetric function, 
restricted to the interval |r| < a, while the window is set 


Gemumity at the origin. 


f £ WCF) df = 1 Co) 


= 00 


w(o) 


w(r) for |r| >a (5.2) 


H 
O 


The frequency response of this window is given as: 


Wer) = r wlr) J_ (fr) dr (5.3) 


Os 


The frequency response consists of a main lobe and side 
lobes which posses a large amplitude ratio in a linear scale. 
The frequency response plot is given in log-frequency scale 


versus log-normalized amplitude (W(f)/W(0)). 
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A. RECTANGULAR WINDOW [2] 
The two dimensional version of the circularly symmetric 
rectangular window is given as follows: 
1 for |r| < a 
wir) = { (5.4) 
0 for |r| > a 
This window is cylindrical, with the height equal to one and 
the diameter of the base equal to 2a. The vertical cross- 
section of this window is equal to a rectangular window. 
The frequency response of this window is given by the 


Hankel transform as: 


a 
W(f)- fr J (fr) dr Sa) 
O 


The recurrence formulae is used to solve the integral 


MOS. da), where: [10] 


É GJ (Q0) 


l 
X 
- 
C4 
X 
— 


+ {x J4 GO] =: J (x) 


Then the integral form in (5.5a) reduces to: 


a 
WE) = H 4 (fr)J, (fr) d(fr) 
Q 
pO o 
a a 
= 3 / d(fr)J (fv) = Le | 
O P O 


66 





J. (af) 
A AL, 
- à ae e (Ss) 


J, (af) 
The window frequency response ED has a form similar to 


E form. The value of this function at the origin, 
f=0, is one. The envelope of the window frequency response 
crosses the frequency axis, whenever J, (af)=0. The first 
zero-crossing occurs at the value of af=3.8. If the value 
of a is increased (meaning the window impulse response is 
lengthened), the zero-crossing occurs at smaller value of 
f (meaning the window frequency response is narrower). If 
the value of a is decreased (meaning the window impulse 
response is shortened), the zero-crossing occurs at larger 
value of f (meaning the window frequency response is wider). 
At the value of a=0.5, the first zero-crossing occurs at 


Re value of f=7.6. 


The ENBW of this window with the value of a is given as 


Bene nl) = 3.14 


The COH.GAIN is given as: 


Ds US 
COH.GAIN » ———4 /| P dp S aa 


SUN s O 


tt 
|H 
e 


The window and its frequency response are given in 
Figure 5.1(a) and (b) respectively. The main lobe width 
(the width between the origin and the first zero-crossing) is 
7.6. The first side lobe level in 18 dB down from the main 


lobe peak and the side lobe fall-off rate is 9 dB/octabe. 


E 
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- Figure 5.1 


(a) 2-D Rectangular Window, (b) 2-D Spectral Window 
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B. BARTLETT WINDOW 
The two dimensional circularly zymetric Bartlett window 


1S given as: 


mor) = 1 (5.7) 
0 elsewhere 


The shape of this window function is similar to a cone, 
with the height equal to one, and the base in a circle with 
a radius r=a. The vertical cross-section of this window is 
equal to a one-dimensional Bartlett window. 

The frequency response of this window is given by the 


Hankel transform as follows: 


a 
WCf) mer] Ta J (fr) dr 
O 


a e 
= f pJ (fr)dp - / EJ (fr) dre 
O a o 
O O 
l) J, (at) 1 4 d (fr)J, (fr) 
deca 7 2 
O de 
J. (af) Urea a (fr)J (fr) 
a,» 
E O f 
2 

1 a (fr) Ji Cfr) 
= 2 f m d(fr) 

a y 

O er 
2 

1 4 d(fr) y (fr) 
Z 5 Of 

a 3 

O frp 
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a 
En 1 2 1 
O f r 
J, Caf) 1 a 
E T a d J,(fr) dE) 
f af O 
J, Cat) 1 a 
di P | i J (fr) alp) = 24, laf)) 
f af O 
er LS 
er: {af J, (af) - 2 J, (af)) N Cor) dd (i 
O 
af aro 
1 E a 
== (-af J_(af)} + 3 fr J Cfr) d(fr) 
af (af) O 
D d = 
ERO) = a (1 x J (af) + 3 a ir) d(fr)} (5.8) 
(af) S (af) o 9 


eai rs substituted into (5.8) then the end result is: 


se =) 
f J (t)dt - p J Co)) (5.9) 


Wp) = = 
O 


where t is a dummy variable of integration. 
The window and its frequency response is given in Figure 
5.2(a) and (b) respectively, where a=0.5. The ENBW and 


COH.GAIN parameters are given as: 


ENBW = n(1.33) = 4.13 


27 Diss 
COH.GAIN - ae AS 
M COS?) O 


O 

© 15 

o1 

ft 
ar 
P 
| 
O 
N 
C» 
IH 
O 
GO 
"m 
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The main lobe width is about 11.5, The first side lobe 
level is 34 dB down from the main lobe level, and the side 
lobe fall off rate is 12 dB/octave. 

The term ` 3 (tat in (5.9) is hardly a closed form, 
But it is Pibutated: Then the rectangular and relatively 


RS Bartlett windows are two-dimensional windows that have 


a closed-form frequency response. 


C. COSINE WINDOW 

The third two-dimensional window to be discussed is a 
cosine window. The two dimensional circularly symmetric 
Window is given as: 


cos r on a 
a — 


Mr) = 1 (5:10) 


O elsewhere 


The shape of this window is similar to a half sphere and 
the height at r=o is one. The base of the window is a circle 
with radius r-a. 

The frequency response of the window is given as a 


Hankel transform with no closed form. 


a 


mor) = J r ceos = J (fr) dr (Sell) 
a a O 


The effort to derive a closed form is given in the Appendix. 
If a=0.5, then the window and its frequency response are 


given in Figure 5.3(a) and (b) respectively. The parameters 


I2 
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Ereune. 5.3 


(a) 2-D Version Of Cosine Window (cross-section), 


(b) 2-D Spectral Window 
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ENBW and COH.GAIN are given as: 


EXNIBW"= T(1.33) = 3,86 


OT DNS e 
COH.GAIN » ————— [| r cos Tr dr = == = 0.46 


(0.5) O 78 


©} © 
Go 


im main lobe width 1s 10.5 while the first side lobe 
ers 26 dB down from the main lobe and the side lobe 


Men off rate is 12 dB/octave. 


D. HANN WINDOW 
The two dimensional circularly symmetric Hann window 


is given as: 


cos“ ze Le er Ip asa 
( a a 


w(r) = (5.12) 


O elsewhere 
The vertical cross-section of the window is equal to a 


one dimensional Hann window. 


The frequency response is given by the Hankel transform 


O 
f r(0.5+0.5 cos = p) Jo Cfr) dr 


ar) = 
a 
J, (af) a E 
Es T n.5 JJ cos —r J (fr) dr Comes) 
iE A a O | 


74 





The first term in (5.13) has a closed form as in the case 
of a Rectangular window. The second term has no closed form 
as in the case of the cosine window. 

For a window length of a-0.5, the window and its frequency 
response are given in Figure 5,4(a) and (b) respectively. 


The ENBW and COH.GAIN parameters are given as: 


ENBW = m(1.5) = 4.7 


eas 
m se asso medr = nn S ees 


"ONE O 


O 


COH.GAIN = 


E main lobe width 1s 13.7 and the first side lobe level 
is 34 dB down from the main lobe rele ori rate 


Gr che side lobes is 18 dB/octave. 


E. HAMMING WINDOW 


The two dimensional circularly symmetric Hamming window 


1s given as: 


0.54 + 0.46 cos =P lr] <a 
wir) = { (5.139 


O elsewhere 


The shape of this window is a smooth cone sitting on a 
cylindrical base. The vertical cross-section of the window 


is equal to the one dimensional case. 


The frequency spectrum of the window is given by the 


Hankel transform as: 


T5 
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Figure 5.4 


(a) Von Hann Window in 2-D, (D) Spectral Window 
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a 
ee) fF Oe. O.o4 + 9.56 cos Ap) J (fr) dr 
O 


Ji Caf) a - 
290-513 n LN be nm eos SE J (fr) dr (SO) 
O 


The first term in (5.15) is the familiar frequency re- 
sponse of the rectangular window and the second term is the 
response of a cosine window. 

If the length of the window is a=0.5, then the window 
and its frequency response are given in Figure 5.5(a) and (b) 


respectively. The ENBW and COH.GAIN parameters are given as: 


ENBW = 7(1.36) = 4.3 


OT Us 


MET] = el ; r(to.54 + 0.46 cos 2) dre 
> DR 
7(0.5) O 


062028 
gs 





-0.36 


ime main lobe width is equal to 13.5 and the first side 
be" level 15 47 dB down from the main lobe level. The side 


lobe fall off rate is 6 dB/octave. 


F. BLACKMAN WINDOW 
The two dimensional circularly symmetric Blackman window 


1s given as: 


DAS cos— + BIB cos r In] «a 


wer) = 1 (5.16) 
O elsewhere 
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(b) Spectral Window 


(a) Hamming Window in 2-D, 
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The vertical cross-section of this window is equal to 
the one-dimensional Blackman window. 


The window frequency response is given by the Hankel 


transform: 
a T 2T 
W(£) = fJ r{0.42+0.5 cos- r * 0.08 cos r}J_(fr)dr 
O 
Ji (Caf) a m 
= Oct 2 a —+0.5 J r cos— J (fr)dr 
E a ^o 
= 2T | 
IAS MA cos D J (fr)dr (5.17) 
O 


If the extent of the window given as a=0.5, the window 
and its frequency response are given in Figure 5.6(a) and (b) 


respectively. The ENBW and COH.GAIN parameters are given as: 


ENBW = m(1.57) = 4.9 


DSS A 
COH.GAIN = NN UA. ADA OU + Ones E MUR NOB cos T d. 
2 Ol 0725 
Geo) o 
| er > 
= 7-75 = 0.22 


The main lobe width is equal 18.5, which is almost triple 
of that of the rectangular window. The first side lobe level 
is 60 dB down from the main lobe and the side lobe fall off 


pate is 9 dB/octave. 
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Figure 5.6 


(a) 2-D Blackman Window, (b) Spectral Window 
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G. KAISER WINDOW 


Two dimensional circularly symmetric Kaiser window is 


given as: 
I [8 A - (125 
O a Ir] <a 
_ ieee = 
ar) = 1 O (5.18) 
O elsewhere 


The vertical cross-section of the window is equal to a 
one dimensional Kaiser window. 
The frequency spectrum of the window ls given by the 


Hankel transform: 


I [8 V1- Reb 


W(f) = E — Jo Cie dae Goo) 


O DU 


which cannot be reduced to a closed form. 
If a=0.5 and 8=6 then the window and its frequency res- 
ponse are given in Figure 5.7(a) and (b). The ENBW and 


COH.GAIN parameter are given as: 


Ben = (1.5) = 4.71 
——— 
T 0.5 I Bir Ez ] l up | 


EN ug. pr n Sp e < 
Deo 15:8 0.78 


© 


The main lobe width is 14,0 and the first side lobe level 
is 44 dB down from the main lobe peak. The side lobe fall off 
meus is 6 dB/octave. 

In this case, the window's characteristics depend strongly 
on the values of a and 8. With 8-6, the window is similar to 
a Hamming window. 
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(a) 2-D Kaiser Window, (b) Spectral Window 
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TABLE VI.I 


Two-Dimensional Window Functions and Figures of Merit 


WINDOW ENBW COH.GAIN a] a, b d 
dB dB H7 dB/octave 
Retangular 3.14 ISO -18 -38 ON -9 
Bartlett 4.13 . 34 -34 -58 1.140 -12 
Cosine 9236 .46 -26 -55 LOO -12 
Hann 4.71 29 -34 -64 jest -18 
Hamming 4.3 23D -47 -51 C -6 
Blackman 4.9 22 7 -60 -66 Ne? -9 
Kaiser 4.7 or -44 -61 14.0 -6 
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VE  CONCEUSTON 


One of the most straightforward applications of the rec- 
tangular window is to truncate an infinite length impulse 
response. This can be thought as multiplying the impulse 
response by a window function that is unity in some finite 
interval between ta, and zero elsewhere. This multiplica- 
tion is equivalent to convolving the Fourier transform of 
the impulse response with a function of the form m This 
causes ringing in the original frequency response. Increasing 
the length of the finite impulse response (by increasing the 
length of the rectangular window) decreases the duration of 
the ringing, but has no effect on the peak overshoot. This 
effect is called the Gibbs phenomenon. 

A considerable amount of effort has been put into finding 
wea window function, i.e., window functions which reduce 
the peak overshoot of the Gibbs oscillation. It is generally 
agreed that the frequency response of a good window should 
have a narrow central lobe (with a large amplitude) and side 
lobes with comparatively small amplitude. There are many 
window functions readily available in the literature, however 
most window functions are given only in one dimension. The 
one dimensional window functions and their Fourier transform 
are given in detail in Chapter III of this thesis. Figures 


of merit for one dimensional window are powerful measure and 


are used to compare the characteristic of the various windows. 
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It has been the objective of this thesis to examine two 
dimensional window functions. The two dimensional window 
functions presented here have been a direct extension of 
corresponding one dimensional window functions. The two 


Gammensaonad window function is formed from a one dimensional 





window function, using the relation uo y eu Oa y^). This 
obviously rotates the one dimensional function about the 
origin in the x-y plane. However it does not simply rotate 
the frequency response of the one dimensional about the origin 
in the u,v plane. The two dimensional frequency response of 
the window can be calculated using the Hankel transform. The 
principal advantage of this approach is that the frequency 
transform of two dimensional window can be found using a one 
dimensional manipulation. 

The two dimensional window figures of merit are a direct 
extension of the one dimensional case. However the definition 
of the ENBW is changed from the width of a rectangular filter 
in the one dimensional case to the circumference of a circular 
filter in the two dimensional case. Table VI.1 shows the 
comparison between one dimensional and two dimensional window 
figure of merit. From the table it is clear that the two 
ENBW parameters are related by the equation (4.42). 

The coherent gain decreases slightly except for the Kaiser 
window where it is higher, and the rectangular window where it 
is constant. The first side lobe level, which is shown by the 


parameter a is lower for the two dimensional windows, however 


ng? 


the width of the main lobe is much greater in each case. The 
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side lobe fall-off rate stays the same except for the Blackman 
window where it is slower and the rectangular window where it 
is faster. 

In conclusion, while one dimensional and two dimensional 
windows can be simply related to one another in the time 
domain, no direct, general parallels exist to relate their 


transforms in the frequency domain. 
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CSE E TL 


Comparison Between Two Dimensional And 
One Dimensional Window Figures of Merit 


ENBW COH.GAIN a] b 
WINDOW 

2D Mp 2D 1D 2D 1D 2D 1D 
Rectangular 3.14 1.0 i. cl MS SS SA 0. 
Bartlett DEN 33 0 U n5 DESDE lI. 
Cosine Secor. cs «tb 0.64 -26 =23 10.01. 
Hann E LS 0 Use COM  S2 e 1: 
Hamming 4.3 eb Faser), Su 2ye L3 28.01. 
Blackman t1 57 22 SIPESESSEDOTU Berl Ek 2. 
Kaiser 4.7 uS qb .49 -44 -43 14.0 2. 


87 





APPENDIX 


The effort to put a Hankel transform of a cosine function 
in a closed form is given below. The cosine function arises 
from Hankel transform of a cosine window, 

The Hankel transform of a cosine function is given below: 

a 


w(t) = f r cos 
O 


T 
zá P JPE] dr Ir! Sl CASS) 
The above integral form has no known closed form, as proved 
by the following discussion. 

Three different methods are used to bring the integral 
Eua closed form, they are: straight integration, partial 


differential and table of integration. 


A. STRAIGHT INTEGRATION 
To use straight integration method to solve the integra- 
tion, the variable needs to be changed into trigonometric 


Ewsbles. 


The interval of the integral becomes, 0-o and =z for r=o 


hoya 


and r=a respectively. 


Ms titution of the new variables into Hankel transform 


gives: 
1.2 7/2 
RAGE) = (7) if 6 cos 5 JotA9) dé (A72) 
O 
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where À = US. This new integral form gives more pleasing 
appearance than the previous one. 
Substitute the Bessel function J¿(A0) inside the integral 


with Bessel equation, given as follows: 


i 


Bey) 4$ —J ' (x) +J (x) = 0 CARS) 
O O O 
Let x = 10, and do some manipulation. The Bessel equation 
becomes: 
Mee) = A e d '(9)) (A.4) 
O Y O Eno S 


The integral form in A.2 becomes: 


MZ 


WA) =-5 J0cos 0 (UJ "(10) + 2 J '(A0))d6 

ils O S = 
1 TZ 2 

= s casa UCNO) | f [(cos 8-8 sin 8)J '(498) 
f z o o S 

- cos 0 J '(à90)}d9} 

1 1/2 

NM sino J-'(X0)- de) 
f O 2 
1 1/2 m2? 

o o sân 6d (19) | = f (sin 6*9 cos 8)J (A6)0d60) 
E E O O = 

1 2 
E TT T l a É J 
- uo DEO - = : S J (A98) de + 1 W(A) 
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at = ZA) = 
\ 


The end result is: 


w(A) 


f sin 6 J i648) do - 


19722 


T TT . 
EE uM UT — y sin 9 J (A8) dé 
ago 9 2 E 2 


2 
Ex J. (GA) (A.5) 
2)“ E E 


By replacing J (AB) inside the integral as before, 1t may 


be continued at will, but the term with the integral form 


cannot be eliminated. The conclusion is, the integral form 


in A.2 cannot be put in a closed form by this method. 


Pee ant IAL DIFFERENTIAL METHOD 


Return to the integral form given in A.2: 


WOA) 


Define 


CLA) 


SCA) 


| > 


1/2 


0 cos 6 J 00) d8 


new equations: C(A) and S(A) 


9 cos 9 J (AB) do (A.6) 


8 sin 0 J (A) do (ATT) 


M ferentiate to A as C'(A) and S'(A) 


C'(2) 


n/ 7 > 
= Ie O J (AB) do (A.8) 
O 
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a 


S' (A) = B 


SCRO m £f S2 ERG J '(A8) de 


Multiply both A.8 and A.9 with A: 


NEN) = f 0% cos 6 ren de) 
Q 
TZ 2 

ASICA) = f o° sin 0 J,'(A0) d(A0) 
O 


(A. 9) 


(A.10) 


GAR) 


Using Bessel identity Jo (A )d(AB)=dT (A9), substitute 


to A.10 and A.11l: 


1/2 


AC! (3) NCC daJ (A0) 


| 
—, 


RS (A) 


l 
— 


o? cos 8 dy (49) 


1/2 1/2 


(RE) 


CAS) 


pana both integral in A.12 and A.13 into partial integration: 


AC'(A) - 8* cos 8 J (A0) = ue (28cos9-0 sine)J (18)d0 
O O O O 
m. 2 2 
A cos0-0 sine) J O48) do (A. 14) 
O 

2 1/2 2 

AS'(A) z 0" sin6 J| (A0) | (20 sin8*6' cos0)J (A0)d60 
O 


SE: 





0/22 
JG - f (28 nen cos0)J, (A0)de ees) 
Q 


i 
ia 


Rewrite both equation in A.14 and A.15, using definition 
EN sand S. It will end up with the following pair of differ- 
ential equations: 


NOT CA) a + SECA) (A.16) 


AS" CA) A) =- 25 - C! (A.17) 


o/a 
C, 


( 
O 


N| > 


Rearrange the last two equations: 


AE" + 20 + st = 0 (A.18) 
mae cos ct = 3 (da) (A.19) 
2 O 2 


The last two pairs of equations show that it may give a single 
equation in terms of S or C. And that means it will end up in 
integral form of C or S. The conclusion is, the integral 


form in A.2 cannot be put in a closed form by this method. 


C. TABLE OF INTEGRATION 
Taylor series expands any analytic function f(x) into a 
power series of x. Consequently cos x may be expanded into 


2 
a power series of x 
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Born = Ce (A.20) 
Substitute 4,20 into A,2: 
f 2 
mex) = J x f(x) Jo (fx)dx (A.21) 
O 


where we change the integration variables with the dummy 
variables x and f. 

From Table of Integration [Ref. 10], the integral form 
shown below: 


¿2K+1 J (tat (A.22) 


O =A 


Thy (x) = 


can be written ina closed form. The term t(3,(t) in the 
above equation is replaced using Bessel equation: 
meet) = -tJ "(t)-J '(t). 
O O O 
The integral form in A.22 becomes: 


= E 2° t t 
ne.) amen tet NE) = I 'CEP)) dt 
2k is O O 


2k+1 SE 


k 
--f t EE CE) ae 
O 2 


xX k 
-x KTIJ 1 (x)+(2k+1) f t^*j '(t)dt- T t? 


O O 


Ky UT Gap reer 
o 
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` x 
-x KLJ t (x)+2k{x “J _(x)-2k fr‘ 


O 


“ly (t)dt) 
O 


NU (x) = 2kx 5J (x)-x ET 1 (o car 


>k (x) (A.23) 


2k-1 


The last equation shows that I 4X) reduces to lower order 


2 


I (x). This form is called reduction form. For k=o, 


2k-l 
A.22 and A.23 become: 


X 
ieee) = JJ Ct)dt = =xJ_'(x) (A.24) 
O t O O 


The equation in A.24 agrees with the property of Bessel 
DEEL on, 

If the reduction form applied repeatedly and written in 
a general form, A.23 may be written as: 
k 


X 
Ji it 
O 


tJ (t)dt s 2k.Poly(2k).J, G *x. Poly(2k).J, ' (x) (A.25) 


Lets return to equation A.20: 


Eos x = e) = 5 É 2k R 
k=0 
Em Statute A.26 into A.25: 
e 2 5 2k i 2k*1 
feet ot J (tJdt = X ax de) pee) bo x SCIT Cx) 
O E k O ` k O 
O SO ¡NO 
(A.27) 
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Apply differentiation to both left and right hand side of A.27: 


CO 


CO CO 
R a H IT TT 
p k O O S k O E k 
leo lem lC ZO 
A a „ek EE S x4K*d 
O k O k 
Replace J "O0, using Bessel equation: 
Mey eee a tC) = J Cx) (A.28) 
Ö O O 
NES =) tt ora AT K b th 
E o x O O T k k 
2k 2k 2k 
t 
+ d. (x) (Za, x + E (2k+1)b, x -Eb, x } 
` = 2k+1 — 2k+1 
S xt e (2k*1)a,. x & e bx } 


ren 
O 


+2kb ws 
S k 


k 


Equate the coeficient of the equal term d. and d. from the 


left hand side and right hand side: 
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CO 

2 it 2k+1 

x fX = : t(2k+1)a, 44 bL) X (A.29) 
u 2 
DES E (a, *2kb, )x (CA 30» 
Therefore: 
f +b 
T LN E ~“k-l 
É = (2k*l)a, 17b, > a,” ESAE EET (A.31) 
ax 

H = a *t2kb, > Dy = -= 2i (A.32) 


lu. last two equations give the solution to the integral 


of the cos x form. Lets expand cos x into power series: 


0-9 d fx t t.x t f.x + 
O E 


Een using relations given in (A.30) and (A.31): 


a = 0 b = Q 
O O 
L. Eo 
E - >. Br = = 7 
` EE SC. E 
B^ — n 2 16 
E. f,-f£ 5/4 i RU 
3 E 3 192 
etc. 
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The conclusion is, we succeeded in solving Hankel trans- 
Berm of Cos x function, but we failed to express it in a 


closed form. 
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